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ABSTRACT 


A  program  was  written  for  the  IBM  360  system  of  the 
Naval  Postgraduate  School  Computer  Facility  to  apply  the 
methods  of  singular  perturbation  theory  to  earth-moon 
trajectories . 

Verification  of  results  against  previous  work  was 
carried  out  and  it  was  found  that  agreement  could  be 
attained  only  for  energy  parameter  values  of  0.707  or  less 
No  solution  for  higher  values  could  be  found. 

The  analysis  of  three-dimensional  orbits  was  then  con- 
ducted within  this  restricted  range  to  show  the  merit  of 
singular  perturbation  theory  as  an  initial  design  tool, 
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I.   INTRODUCTION 

With  the  advent  of  space  travel  and  of  complex  satellite 
systems,  a  revival  of  celestial  mechanics  has  appeared  in 
modern  science  and  engineering.   In  its  present  form, 
commonly  referred  to  as  orbital  mechanics  or  astrodynamics , 
it  combines  the  theories  of  astronomers  and  scientists  such 
as  Kepler  and  Newton  with  modern  ideas  and  technology  so 
that  man  may  navigate  his  way  through  space. 

Orbital  mechanics  problems  grow  in  complexity  as  the 

number  of  bodies  under  consideration  grow.   In  fact,  they 

become  so  complex  that  only  the  problem  of  two-bodies  may 

be  truly  solved.   Solutions  to  forms  of  the  three-body 

problem  are  accomplished  under  certain  special  conditions, 

but  usually  numerical  approximations  exist  for  these  higher 

order  problems.   Szebehely  [Ref.  1]  states, 

"Not  only  is  the  two-body  problem  solved  —  and 
the  meaning  of  'solution1  may  be  different  for 
astronomers,  engineers,  and  mathematicians  ■ —  but  a 
general  understanding  exists  regarding  this  dynamical 
system.   The  problem  of  three-bodies  on  the  other 
hand  is  neither  solved  nor  is  the  behavior  of  the 
dynamical  system  completely  understood. 

In  order  to  cope  with  the  problem  of  space  navigation, 
many  methods  to  determine  the  motion  of  a  vehicle  in  space 
have  been  attempted.   The  use  of  two-body  solutions  has  been 
used  in  the  method  of  patched  conic  sections.   Many  pertur- 
bation techniques  have  been  applied  for  two,  three,  four, 
and  more  bodies.   The  most  common  method  of  calculating 


these  trajectories  is  probably  the  use  of  the  digital 
computer.   The  computer  is  used  to  solve  the  initial-or 
final-value  differential  equations  of  motion  in  numerical 
form  and  is  often  time  consuming  and  expensive. 

Due  to  the  difficulty  in  finding  solutions  to  trajectory 
problems,  Lagerstrom  and  Kevorkian  in  References  2  and  3 
formulated  a  method  in  which  they  attempt  to  treat  the 
restricted  three-body  problem  with  a  technique  similar  to 
that  used  in  certain  fields  of  fluid  mechanics.   This  is 
the  method  of  solution  of  singular  perturbations  by  inner 
and  outer  expansions  and  represents  an  attempt  to  obtain  a 
uniformly  valid  analytic  expression  for  trajectory  motion 
in  the  restricted  three-body  problem.   The  work  of  Lager- 
strom and  Kevorkian  is  aimed  at  giving  a  tool  to  more  easily 
predict  the  effect  of  initial  conditions  on  the  motion  of 
a  vehicle  as  it  becomes  influenced  by  a  third  body  in  its 
field  of  space.   This  would  allow  the  determination  of 
"ballpark"  data  for  use  in  more  detailed  solution  programs 
at  a  considerable  savings  of  time  and  expense. 

In  this  paper,  an  attempt  to  use  the  singular  pertur- 
bation method  of  Lagerstrom  and  Kevorkian  to  study  effects 
of  initial  conditions  on  lunar  trajectories  is  made.   A 
brief  discussion  of  Lagerstrom  and  Kevorkian's  theory  is 
followed  by  some  numerical  applications  of  this  theory  and 
the  results  are  discussed. 
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II.   THEORY 

A.   EXAMPLE  INNER  AND  OUTER  EXPANSION  MATCHING  PROBLEM 

Before  formulating  the  singular  perturbation  solution  to 
the  restricted  three-body  problem  it  is  useful  to  look  at 
a  simple  example  of  matched  inner  and  outer  asymptotic 
expansions.   The  example  is  basically  that  of  Prandtl  and 
has  been  used  by  Van  Dyke  [k] ,  Ames  [5]  and  others  as  a 
classical  singular  perturbation  solution. 
If  an  equation  of  the  form 

is  written  and  a  regular  expansion  is  attempted,  difficulty 
is  expected  since  the  equation  obtained  when  M  =  0  is  of 
lower  order  than  the  original. 

Consider  the  following  equation 

/**"■•  r  .  *»  (i) 

with  initial  conditions, 

?l<*  --  o  and        U^  ■=   \ 

The  exact  solution  to  this  equation  is 

Kx.^rO-^l'-e'^J/tl-e'7^  +tx        (2) 

If  a  regular  perturbation  expansion  is  attempted  one  may 
write 
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r.*:*/4>A*M:-  •  •  • 

r--T." - ^ v: - /* U'  ,  ••    • 

substituting  into  1 

matching  powers  of  ^  implies 

t:  - 1 

so  that  by  integration 

$.  =  1*  4  C  (3) 

The  boundary  conditions  cannot  be  satisfied  unless  b  =  1. 
If  the  boundary  condition  at  x  =  0  is  dropped 

since 

^OW  I  =  Uo-vc  — *  C=  i-b 

therefore 

tU;H  =  bx  +  Ll-b^  (i|) 

This  approximates  the  solution  quite  well  except  in 
the  region  adjacent  x  =  0  (commonly  referred  to  as  the 
boundary  layer),  see  Figure  1. 
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Ames  states  that  if  terms  are  lost  or  boundary  con- 
ditions discarded  in  the  outer  solution  they  must  be 
included  in  an  inner  expansion  which  must  be  developed 
(i.e.,  this  implies  matching  in  the  overlap  region). 

The  asymptotic  matching  principle  states  that  the  m-term 
inner  expansion  of  (the  n-term  outer  expansion)  =  the  n-term 
outer  expansion  of  (the  m-term  inner  expansion).   See  Figure 
2  for  a  diagram  of  the  calculation  sequence  applying  this 
principle. 

Returning  to  the  present  example,  the  first  term  is 


f0  =  O-M  +  b% 


for  the  outer  expansion.   The  first  term  of  the  inner, 
expansion  is  found  by  a  "coordinate  stretching."   The 
dependent  and  the  independent  variables  may  be  stretched 
together  or  separately  using  physical  insight  or  trial  and 
error,  however,  the  stretching  must  include  essential 
elements  omitted  in  the  outer  expansion  procedures  and  must 
be  matchable  to  the  outer  expansion. 

In  the  present  example  one  may  desire  to  stretch  only 
the  independent  variable  such  that, 


where 


{c*,>o=  Hx,>->) 


X   =    Vsc^ 
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Using  this  in  the  initial  equation 

<**"         /*    u>c  /-  (5) 

with 

FM  =  0   ,  FOls)  =r  | 
where 

5(^)o(X  -  d  ^ 

The  highest  derivative  was  lost  in  the  outer  problem 

2     2 
therefore  d  F/dX   must  be  kept  for  the  inner  problem.   This 

implies 

SC}A/j*  &  CO  as    ^  -^  o 

but 

s</*)//*  "  °      as   ^  — *  o 

implies  that  the  solution  satisfying  the  inner  boundary  con- 
dition is  a  multiple  of  X  and  cannot  be  matched.   Therefore, 
the  only  possibility  is  that  SC)A  J/**       =  a  constant  as 
/*•— ♦  0,  and  this  constant  is  assumed  to  be  1. 
This  gives 

55?- A  ir  =  k^  (6> 

with 

FCo)  =  0  ,  FO/^  ^  I  ,  ^^  »  /a- 

In  the  inner  variables  X  =  x/^  and  F,  the  following  expan- 
sions are  produced: 

p  =  F.  -t  ^F.  -f  ^  F-u  +  •  •  • 
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*IX    AX   '  AX    7   AX 


5v*        A*1  <*X     r    A>r 


which  when  substituted  into  5  yield 

fHfl  +  i^  =  °  and   £  Co)  -  o 

and  finally, 

F.t^r  A(l-  e"X)  (7) 

where  the  outer  boundary  condition  has  been  discarded. 
To  carry  out  the  matching, 

1-term  outer  expansion  =  (.1-1)4  o  "*• 

in  inner  variables  =  O- M  *  b>*-X 

expanded  for  small  JLa  =  C\  —  t)  +  K^"^ 

1  term  inner  expansion  =  O-b) 

1-term  inner  expansion  =  A  L  I  -  £   ) 

in  outer  variables  =  A  L  \  -   c"  w 

expanded  for  small  v*-  =  A 

1  term  outer  expansion  =  ^ 

matching  first  term  outer  to  first  term  inner  gives 

A--  i-b 

therefore 

ro^O-b)Ci-e"x  ) 

so  that  now 
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to  find  the  second  term  inner  expansion 


therefore 

F-c\-t)ci-er*)  +  /-l^te-x+o  +tx]+  ••-        (8) 

Since  the  exact  solution  for  the  outer  variable  is 

known,  the  matching  from  1  to  2  may  be  skipped  as  it  simply 
yields, 

|-t   ^  I-  b 
Going  now  from  m  =  2  and  n  =  2 , 

2-term  outer  expansion  =   O-lO  +■  h  Y- 

rewritten  in  inner  variables  =   £l-b)+  ^^X 

expanded  for  small    ^  =   Ci  -  W)  +  >^bK 

2  term  inner  expansion  =   Cl-t)  +^bX 

2-term  inner  expansion  =  O-^G-e'^H^^e^-O+^X] 

in  outer  variables  =  0-OG-e*)+^cfe"*Li)UJ:] 

expanded  for  small  y<  -  (\-\>)  — ^C+-tX  +  •  •  • 

2  term  outer  expansion  =  0-(>)  +  l,X  -J^C 

rewrite  in  inner  variables  =  (j-bji-j+bx  -J*-^- 
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matching  these  gives 

bx  =-  bx  -  C 
C  -   o 
therefore 

\:  =  Ci-fcte\-e.~*l  +  >x\o>C  (9) 

This  form  of  solution  is  necessary  when  no  single 
asymptotic  expansion  is  uniformly  valid  throughout  the  field 
of  interest  in  a  problem.   This  problem  may  arise  through 
factors  such  as  two  characteristic  lengths  as  those  of  the 
boundary  layer  thickness  near  the  body  compared  to  the  chord 
length  of  the  body. 

This  long  sample  solution  is  presented  above  so  that  as 
the  solution  of  the  earth-moon  problem  is  discussed  an 
understanding  of  the  procedures  involved  can  be  felt  even 
though  the  actual  solution  is  too  lengthy  to  be  presented 
in  detail.   The  procedures  are  the  same  in  the  orbit  problem 
as  in  the  example  just  discussed  but  more  complicated  to 
carry  out . 

B.   NATURE  OF  THE  PROBLEM 

With  the  method  of  inner  and  outer  expansions  previewed, 
attention  is  now  turned  to  the  trajectory  problem  in  earth- 
moon  space.   It  is  beneficial  to  observe  first  the  general 
problem  which  must  be  studied,  the  three-body  problem. 

The  three-body  problem  is  one  aspect  of  the  n-body 
problem  of  mechanics.   In  earth-moon  space  the  problem  may 
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be  depicted  as  in  Figure  3.   This  figure  represents  the 
restricted  three-body  problem  which  is  the  second  simplest 
of  the  many-body  problems.   The  restricted  three-body  prob- 
lem, however,  is  still  unsolvable  in  a  closed  analytic 
manner.   The  only  n-body  problem  which  has  been  solved  is 
the  two-body  problem.   It  is  therefore  of  interest  to 
observe  the  nature  of  the  three-body  problem  and  determine 
what  information  one  may  find  from  determinable  quantities. 

Using  the  system  of  three  masses  depicted  in  Figure  4 
it  is  possible  to  define  three  position  vectors  in  space 
by  r, ,  ?p ,  and  r~.   Now,  if  the  conventional  methods  of 
vector  mechanics  are  applied  to  the  system  and  the  necessary 
algebraic  manipulations  are  carried  out  it  is  possible  to 
determine  ten  of  the  necessary  eighteen  integrals  of  the 
total  solution  to  the  three-body  problem. 

Writing  the  equations  of  motion  with  reference  to 
Figure  4 


yv\trt    =  <h  rvi,*w  .r,x      +   <h  **.  **-v 


r„ 


r3 

•it 


wui\   -  ^j  wi.»m-uIii   +.   c-j  h^  w*t  lis.  (10) 


Hi  rt\ 


using  the  fact  that 


and     <"i{  «  Tjc 
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and  adding  the  three  equations 


which  yields 

2.YVS.5-    •=  O  (11) 

which  is  the  conservation  of  linear  momentum. 
This  vector  equation  integrated  twice  then  gives  six  of  the 
ten  available  integrals. 
Now  using  the  relation 

on  equations  10  the  following  equations  are  obtained: 


If  the  vector  relationship 

Yc  Ar.    =  -  f j  x.  ?c 

is  applied  and  the  equations  are  added  the  following  equation 
is  found: 
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J 


which  reduces  to 
noting  that 

the  equation 

21  t*U  £*\  *ft  )     —   ° 

may  be  rewritten  as 

which  finally  yields 

2  wavL^i  kxtl]  =    Cct  (12) 

which  is  conservation  of  angular  momentum. 

This  allows  the  determination  of  three  more  of  the  necessary 

integrals  for  the  three-body  problem. 

The  following  series  of  algebraic  equations  will  lead 
to  the  determination  of  the  final  integral  for  the  three- 
body  problem.   Writing  Newton's  law  and  operating  on  it 
one  obtains: 

jit 


l-  mV  ^  =  Jfa  .  dfx 

Using  this  operation  on  Equation  10 
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J    Al  J 


<\i  •  Jirx 


1   wixVvZ  r   C-,  m.v*,     ^  I  AV1  4-  GiWUM-V"  *A 


J  A?  J       ^ 

the  equations  are  now  added  to  find 


r  -  CiVW; 


so  that 


*"M 


(13) 


This  is  the  tenth  integral  which  represents  the  total 
conservation  of  energy.   With  eighteen  integrals  necessary 
for  solution  and  only  ten  integrals  that  are  defineable 
further  pursuit  of  an  analytic  solution  of  the  three-body 
problem  is  impossible. 

If  a  restricted  three-body  system  is  assumed  one  further 
integral  may  be  found  to  offer  new  information.   This 
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integral,  called  the  Jacobl  integral,  does  not  solve  the 
problem  in  analytic  form  however. 

The  following  two  assumptions  are  adopted  to  formulate 
a  restricted  three-body  problem: 

1)  m,  «  m,  and  iru 

2)  m,  and  m^  are  in  circular  orbits  about  each  other. 

With  this  preliminary  introduction  to  the  three-body 
and  the  restricted  three-body  problems  concluded  attention 
will  now  be  turned  to  a  singular  perturbation  solution  to 
the  restricted  three-body  problem  formulated  by  Lagerstrom 
and  Kevorkian  [2] .   The  development  will  be  outlined  in 
notation  of  Reference  2,  however  for  more  detailed  dis- 
cussion and  mathematical  details  the  reader  is  referred 
to  Ref.  2.  . 

C.   FORMULATION  OF  THE  RESTRICTED  THREE-BODY  PROBLEM 

With  a  brief  discussion  of  the  general  three-body 
problem  completed  and  with  an  example  problem  of  matched 
inner  and  outer  expansions  discussed,  attention  is  now 
turned  toward  a  discussion  of  Lagerstrom  and  Kevorkian's 
work  in  the  application  of  this  expansion  solution  to  a 
restricted  three-body  problem  in  earth-moon  space. 

The  planar  restricted  three-body  problem  may  be  formu- 
lated considering  the  notation  used  in  Figure  5.   In  this 
figure  x-y  is  an  earth  centered  frame  of  reference  and  x-y 
is  a  moon  centered  frame  of  reference. 
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Motion  in  the  x-y  system  may  be  considered  to  be 
Keplerian  about  the  earth  with  only  small  perturbation 
terms  due  to  the  moon  until  the  particle  comes  into  the 
influence  of  the  moon's  gravity  field.   At  this  time  the 
perturbation  from  the  moon  is  no  longer  small,  hence  the 
second  characteristic  length  which  results  in  the  singular 
perturbation  problem  and  allows  the  use  of  the  boundary 
layer  "stretching"  type  solutions.   The  effect  of  the  moon 
destroys  the  uniform  validity  of  the  Keplerian  solution 
about  the  earth.   If  the  mass  fraction  of  the  earth-moon 
system  is  used  as  a  scale  factor  one  may  discuss  the  motion 
near  the  moon  in  terms  of  the  inner  or  x-y  stretched 
coordinate  system. 

Considering  the  outer  problem  first  it  is  possible  to 
write  the  planar  equations  of  motion  in  x-y  space  directly 
as : 

(14) 

where 

If  the  terms  on  the  right  hand  side  of  the  previous 
equations  are  set  to  zero,  one  has  the  Keplerian  or  two-body 
equations  near  the  earth.   The  Keplerian  integrals  will 
change  little  near  the  earth.   These  integrals  are  those 
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of  energy,  angular  momentum,  and  position.   Recalling  the 
equations  for  the  three-body  problem  discussed  earlier  one 
has , 

and 

Rewriting  these  in  two-dimensional  coordinates  in  the 
notation  of  Lagerstrom  one  has 

JL[k"+^]  -     »ZT    .ecu    =   Ke  (15) 

and 

X^-^%    =    Co~»7rJe  (16) 

Two  equations  specifying  the  position  of  the  orbit  in 
the  x-y  plane  may  also  be  written  as 

(17) 

These  are  written  such  that  if  p   =  0  the  motion  is 

symmetric  about  the  x-axis.   Also  if  qe >  0  the  x-coordinate 

of  the  perigee  is  negative  for  I     i-   0  and  always  non-negative 

for  £   =  0. 
e 

If  the  initial  conditions  are  specified  in  terms  of 
these  integrals,  (which  will  be  done  shortly)  the  motion  of 
the  particle  in  the  Keplerian  orbit  may  be  determined. 

The  launch  is  considered  to  be  at  perigee,  hence 
p  =  0  is  given. 
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The  angular  momentum  is  defined  also  to  be  rv  cos  <J> . 
If  perigee  is  considered  <J>  =  0  and  cos  <f>  =  1 ,  therefore, 
£   =  rv.   If  this  is  non-dimensionalized  with  respect  to 
the  angular  momentum  of  the  earth-moon  system  (moon's 
period)  the  non-dimensional  angular  momentum  may  be  defined 
as , 

v      rv 

A"  \>!>6,<>e  +  VvO>^  (l8) 

The  non-dimensional  energy  is  defined  in  terms  of  a 
non-dimensional  potential  energy  and  a  non-dimensional 
kinetic  energy  as, 

f  =  \CJ-K>7  -  ^Xx>/Zc^m^m^\  (19) 

Since  h  ,  J^  ,  p   determine  q   the  particle  initial  con- 
ditions are  specified  in  terms  of  h  ,  Pe*jle>  and  q  . 
Writing  these  in  their  non-dimensional  form, 

heV/\    J(e=^"xA  ;    fc  =  O  ^  ^  >  o  and    t=  oC)3} 

one  is  now  able  to  specify  the  initial  conditions  of 
velocity,  altitude,  and  position  in  terms  of  initial 
Keplerian  integrals.   The  following  Keplerian  relationships 
at  perigee  are  known: 

—  xP=  rP  =  olCv  -O 
x  =  o 

C^  =  ftM(Su  -t^ 

Replacing  he  by  -p  ,  £e  by  y 2X  and  pe  by  0  at  x  =  0  and 
carrying  out  some  algebraic  manipulations  the  initial 
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values  for  the  coordinates  and  velocities  become: 
X- 


=  o 


(20) 


where 


and 


4JL 

-   0 

At 

a  - 

0-/-)  /-zp 

-  [i-[w*vo-^1\ 


It  is  also  beneficial  to  transform  the  motion  equations 
to  the  independent  variable  x  rather  than  t.   One  trans- 
formation is  shown  below  and  the  rest  follow  similar  steps 

Writing 


to  produ 

:ce  the 

transformed 

equations 

S-fc 

i  «o  - 

oiy  oU  "Ax 

- 

jt*L  tn  t7 

finally 

= 

HO1    * 

T 

2    2 
A  similar  development  may  be  used  to  transform  d  x/dt  ,  h  , 

and  I      so  that  these  equations  are  finally  available  in  the 

form: 
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D.   OUTER  EXPANSION 

Attempting  a  solution  of  the  equations  of  motion  for  the 
particle  in  earth-moon  space  in  a  power  series  in  k  one  has; 

^lX,)*}  ^/^'^utx^+j^vCxU  •  •  •  (22) 

tO,^  ~  -to  L*j>*-)  +  ^t»6c)  4-  -  -   •  (23) 

The  terms  i>^  2y,  and  t   are  the  exact  solutions  for  the 

equations  of  motion  with  the  perturbation  terms  on  the  right 

hand  side  set  to  zero.   In  common  with  celestial  mechanics 

notation,  these  terms  are  written  as, 

i0  z=  G^LfcUO  -e  s„otCx\|  (24) 

/'XVK.-     -t  S„o  ECO  (25) 


where 


EC*)  =  s^-'Li-Ce-fcyn1* 
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If  toCx,/*)   and  y<   ^.btyOare  further  expanded  as 

£oC*jy)    -  too  C*)  +■  ^V0»  LO  +  •  •  •  ^26) 


^O,^  =  ^i,J.*)   *  ^^-VJ^Cx^  •  •  •  (27) 

where  by  appropriate  expansions  in  powers  of  U  one  finds, 
tooCx)  ^O^pCx-  fWv-  ^pCx-r^V'x]  /  ^Y3         (28) 

£o,^^-t.^)+tLc\-^^xx+iKvyzc^xei-^x^"'u       (29) 

li/zobOr  Z^Hx-^x^''1  (30) 

^«,Z|  oo  -  -*pC?_*  ^Xt-z^/zf  ^-/'V'-yi''7-     (3D 

It  is  also  necessary  to  insure  that  the  moon  and  the 
particle  pass  within  the  scale  distance  tA.  of  each  other. 
The  phase  constant  expansion  which  includes  this  condition 
in  the  analysis  is, 

t=t;  +  ^,/lTwx  +^~r;  (32) 

where 

To  -  £oo  C0 

~T*t*  -   -V  ^ 

and  T-  is  an  arbitrary  term  which  assures  that  the  moon  may 
be  allowed  to  interact  with  the  particle  as  desired. 

The  perturbation  terms  t-,(x)  and  y-^x)  may  also  be 
determined,  the  details  of  which  are  long  and  tedious  and 
may  be  found  on  pages  689-690  of  Shi  and  Eckstein  [6].   The 
final  form  of  these  equations  is : 


-Loo  [etLX)  -  *Liki\  A^fc 

(33) 

^,)  =  -[^0-^"Vz  +l^^^-^)*^T»,ll/tV       (35) 
also   f    (z)    and   gQ(z)    are   given   in   terms    of  the  moon's   posi- 


where 


and 


tion    coordinates    as, 

(37) 


tw*. — -r1 -ns,  -  Usb^  -T.l 


j^  =    s^L-b-co-T,]-^.  _  „„  [^^      (3B) 

If  t,(x)  and  y-,(x)  are  now  expanded  asymptotically  as 
x  -*■   1  the  following  expansions  are  produced: 

■L,U)^    L?>-T-f>'l-)~XAo<tL\-x)   -  Y(^  loco  (39) 


«J.tiO  =  C3^plf^^G-^  f  S(frt  +oO)  (HO) 


where 


r-Jjo'^u)Lo<co-^c^+  5^^  *" 


(Hi) 
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S-\l^WoC^>c^Lpco^co^^%^  4*  4-1      {k2 


with 


Q   =  £j  +V'*'?7-  (i|3) 


and 

u-  i/ie'.co  »*) 

With  the  terms  t  ,  t,,  y,,  and   y  2y   now  available  the 
outer  expansion  for  t  and  y  is  known  and  is  ready  to  be 
matched  near  the  moon. 

E.   INNER  EXPANSION 

As  the  vehicle  approaches  the  moon  it  comes  into  a  hyper- 
bolic encounter  with  this  body.   The  hyperbolic  encounter  may 
be  expressed  in  terms  of  inner  variables.   These  inner 
variables  are  written  as,  .. 

><  =  *~  io*6t  -~)  -  _  y  Si,o6fc-T  )      1   -  t  -To  -y-  ? 

where  y.  is  the  scale  factor  for  the  inner  variables  and  t 
is  a  constant  which  must  be  determined. 

The  inner  variable  equations  of  motion  are  written  as: 


-C5  +§-3  -  -  M^t^l  S/^I  ^a-^.-f    (J|6) 
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where 


r*  =-xS  ^ 


In  the  hyperbolic  encounter  one  may  express  the  Keplerian 
solution  in  parametric  form  as, 

£  =  gL^Cc  s/*A^  -lO)  (49) 

The  Keplerian  integrals  previously  discussed  expressed 
for  inner  variables  are 


K*tC^*^l- 


(50) 

f  =  z£  +  ty? 

These  integrals  are  related  to  the  hyperbolic  elements  by, 


jf  =  -E  s>~B,  % 


e  c©^>5 


To  carry  out  the  matching  of  the  inner  and  outer  expan- 
sions the  values  of  y  and  £  as  functions  of  x  along  the 
approach  asymptote  must  be  determined.   These  functional 
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values  may  be  determined  by  allowing  u  -*  °°.   As  u  -*  °°  the 
following  functional  forms  of  y  and  t  are  determined  to  be: 

l  =  V4  a^u^r-  +  0(a'°  (52) 

The  U1  and  V^  terms  are  the  velocity  components  on  the 
approach  asymptote. 


where 

and  A  and  B.  are, 

A*   -  f/^-t         and    Is  -p/2.tv 

If  the  negative  y-intercept  is  considered  to  be  K,  , 
that  is  to  say  y  =  -  K,  for  x  =  0, 

K,  -  AjA-jjJJ.  (53) 

and 

jf  =■  K,  0,  (54) 

U, ,  and  V,,,  the  velocity  components  along  the  departure 
hyperbola,  may  be  written  also 
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F.   MATCHING  OF  THE  INNER  AND  OUTER  SOLUTIONS 

The  matching  of  the  inner  and  outer  expansions  must  be 
carried  out  at  this  time.   Returning  to  the  example  dis- 
cussed, the  principle  behind  this  matching  may  be  recalled. 
The  only  difference  in  the  two  problems  is  that  in  the 
present  case,  since  only  the  first  term  of  the  moon  inner 
variable  expansion  is  needed,  it  is  only  necessary  to 
rewrite  the  outer  expansion  in  inner  variables  and  then 
carry  out  the  matching. 

For  the  particle  near  the  moon  the  x-coordinate  is  £  . 

This  occurs  when  t  =  T  .   If  the  x-coordinate  is  specified 

o 

to  be 

r.=  l-  J±T,^ 

this  implies  that 

f~=  l-^T.L 

The  phase  relation  was  written  as, 

T=  To  -t  j^Tiu*  >-T,  -Tot^V-i^       (55) 

but  with  t  =  TQ  the  equation  may  be  written  as, 

so  that 

•y  =.  %    +  oco 

with 
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therefore 

r 

If  the  outer'expansion  is  written  as 

t  =  £c»  6)0 +-/*.£  to  I  ^0  +-fclCjc^  +  •-•  (^ 

^  =  /^  ti/J*>°^  *  ^V^  *  -  •  '  (57) 

and  replacing 

<y,    \>^,   1L  -  1  -t~  ^  C  X  -  T»/\  j  -O 

+  Y(p^>]  +■  °c^ 
t«» to  =  To    ***   *-  ^~To  ~/^ 

also  Jut  Is  recalled  that  52  was  written 

l=  i^+ ;.%, -f£-£  ^--o 

The  term  (1  +  U2)"3' 2/°<^  does  not  occur  in  the  later  equa- 
tion.  It  must  therefore  be  cancelled  by  some  other  term. 
Since  the  only  arbitrary  term  in  the  inner  expansion  for  t 
is  ?T   >  a  term  in  t   must  cancel  the  log  p.  term,  hence, 

t   ~   Cl  +  t>^TV"  locf^     +-  C,  (59) 


and 


£  =CT'X  f-O+U^^t-x)  +-b,,Cft--T^  +  Y10-Z;  4-oC,^ 


0) 
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y   may    likewise   be    determined   to   be, 

i_ 

n    ^    }br    -  -to, CO  +~n   +  Ur^-YtM  -*    +oc0 

Comparison  of  51  to  61  and  52  to  60  shows 


If  the  equations  are  written  such  that, 


where 


where 


L=to,CO  -  Hl2r      +   Yip)     -   C\ 


and 


t^e 


then 


(61) 


.ons  ; 

c7,x  ~ 


and 

M  =  -K  -  g-  +  ocn  (63) 


-K^T^  -4o.CO^T  +  S(p>  -YCp^  (61|) 


H»  -  CM   -Bo»,^/(J,  (65) 

Then  it  is  evident  that 

K,  ^K 

Rewriting  t  as 


(66) 


£  -  T£-^L  ^  a3/t/o yf^-.  +  °^')       (67) 
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k+*^l~n>; 


v-,_  tr 


which  completes  the  equation  matching  process  and  determines 

the  unknown  constants. 

These  values  are  used  in  equations  50  to  obtain  from 

substitution 

Ztv  =  I  -+  U  "*■ 

I  -  KU 


\    =  KU  -  u  O  +  uV^ 


(68) 


The  term  t,  may  also  be  completely  determined  to  be 


tr(  -t'cM-I^  fV(P^  +1  -^^gvTfJF"      (69) 


The  total  or  composite  expansion  solution  is  the  outer 
expansion  plus  the  inner  expansion  minus  the  common  part. 
It  may  be  written  as; 

^  *  /^vJv^+^V'O  +^tcyttWoL*«>'l  +  o^  (70) 


(71) 


where 

**  =  ^7  2  -  £  eCi/,co5&  _0,  *.^5  }/uf  (?2) 

f>*  =  (x  -  a  e  Co^  a  }/(J, 

.The  solution  to  first  order  after  moon  passage  may  be 
obtained  by  a  similar  matching  procedure  as  that  used  above. 

Before  leaving  the  analytic  theory  it  is  wise  to  note 
that  this  procedure  of  inner  and  outer  expansions  may  be 
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carried  out  to  higher  orders  of  y.  which  increases  the 
accuracy  as  the  value  of  the  parameter  u  increases. 
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III.   THE  COMPUTER  PROGRAM 

The  computer  program  used  in  this  paper  was  written  for 
the  IBM  360  system  at  the  Naval  Postgraduate  School  Computer 
Facility.   It  has  been  run  successfully  on  the  OS/CMS  time 
share  remote  terminal  and  the  standard  batch  processor  unit. 

The  program  is  a  straight  forward  computation  with  the 
use  of  Scientific  Subroutine  Package  programs  from  the  sub- 
program library.   These  programs  are  used  to  evaluate  the 
integrals  which  appear  in  the  solution  near  the  moon. 

All  computation  is  carried  out  in  double-precision 
notation  to  maintain  as  much  accuracy  as  possible.   The 
integration  routine  is  a  double-precision  sixteen  point 
Gaussian  quadrature  from  the  SSP  library  in  the  Fortran  IV 
language . 

The  necessary  inputs  to  the  program  are  the  non- 
dimensional  boundary  conditions  discussed  in  this  paper. 
From  these  parameters  the  program  calculates  the  necessary 
constants,  then  the  integrals  of  the  motion  equations  are 
computed  for  x  =  1  (near  the  moon).   With  these  values  the 
program  evaluates  the  motion  In  the  vicinity  of  the  moon 
and  prints  the  necessary  parameters.   Control  is  then 
returned  to  the  input  mode  and  a  new  series  of  calculations 
may  be  carried  out.   If  p  becomes  too  large  the  machine 
proceeds  to  stop  and  ends  processing. 
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The  necessary  equations,  a  test  case,  and  the  computer 
program  are  contained  in  the  Appendix.   A  flow  chart  may 
be  found  in  Figures  17  and  18.   All  equations  are  explained 
in  References  7,  8,  and  9  and  the  same  notation  is  used  as 
much  as  possible. 
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IV.   DISCUSSION  OF  RESULTS 

A.  VERIFICATION  OF  COMPUTER  RESULTS  AGAINST  THOSE  OF 
REFERENCE  7 

A  computer  program  to  solve  the  restricted  earth-moon 
three-body  problem  was  written.   The  first  test  of  this 
program  was  against  the  data  presented  in  Reference  7. 
Figures  6,  7,  and  8  show  the  agreement  of  the  present 
results  with  those  of  Ref.  7.   The  results  for  y   and  6  also 
agree  with  Table  1  of  Reference  7.   As  a  final  check 
examples  one  through  nine  of  Reference  7  were  tested  and 
comparison  of  results  was  good. 

No  data  of  any  form  for  values  of  p  greater  than  0.707 
could  be  made  to  agree  with  the  published  data  of  Reference 
7.   The  inability  of  the  program  to  predict  accurate  results 
for  this  case  of  p  greater  than  0.707  was  unresolved.   The 
probable  cause  of  the  instability  of  the  solution  is  due 
to  the  denominator  of  the  fn(x)  term  rapidly  vanishing  near 
this  point.   The  equation  for  fQ(x)  was  rewritten  in 
several  forms,  all  to  no  avail.   The  equation  was  also 
expanded  in  a  Taylor's  series  about  0.   This  expansion  gave 
some  improvement  in  the  prediction  of  y   and  6,  but  the 
improvement  was  not  sufficient  to  be  used  for  numerical 
work. 

B.  PARAMETRIC  STUDY  IN  THREE  DIMENSIONAL  EARTH-MOON 
TRAJECTORIES 

The  program  used  to  predict  the  results  of  Reference  7 

was  modified  to  include  three  dimensional  orbit  consideration 
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for  values  of  i  =  0(1).   The  modification  still  led  to  the 
same  results  for  planar  orbits  and  a  study  for  inclined 
orbits  was  carried  out.   Cases  were  run  to  predict  the 
effect  of  variations  of  energy,  angular  momentum,  position 
of  the  moon  and  orbital  inclination  on  the  space  vehicle 
as  it  approached  the  moon's  position. 

The  present  program  predicts  all  of  the  important 
parameters  near  the  moon.   The  effect  on  changes  in  initial 
conditions  can  be  observed  in  the  changes  in  these  near 
moon  parameters.   One  of  the  most  important  of  these 
parameters  and  one  that  requires  knowledge  of  the  other 
parameters  is  the  parameter  of  closest  approach  to  the 
moon.   This  parameter  has  been  chosen  to  show  the  effects 
of  varying  initial  conditions  on  the  moon  encounter. 

Figures  9  and  10  show  how  the  initial  launch  conditions 
are  related  to  altitude  and  velocity  at  launch  from  perigee. 
The  altitude  limit  on  perigee  shown  in  Figure  9  indicates 
values  of  lambda  must  be  less  than  3.0.   Figure  10  shows 
the  velocity  at  perigee.   Since  conventional  launch  vehicles 
cannot  attain  values  of  velocity  much  greater  than  around 
lJ0,000  kms/hr  an  outside  limit  on  velocity  was  taken  to  be 
50,000  kms/hr.   This,  then  implies  that  lambda  should  be 
greater  than  1.0.   For  a  given  launch  condition  the  values 
of  rho  and  lambda  can  be  related  to  the  velocity  and  alti- 
tude at  perigee.   The  perigee  launch  gives  the  angle  of 
perigee  to  be  JT    and  time  of  launch  to  be  t  =  0.   Hence, 
this  defines  the  launch  conditions  completely. 
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Figure  11  shows  the  effect  of  changes  in  the  initial 
energy,  rho,  relative  the  earth  on  the  distance  of  closest 
approach.   It  is  apparent  that  as  the  energy  is  reduced 
(rho  becomes  greater)  the  distance  of  closest  approach  to 
the  moon  is  increased.   The  effect  of  orbit  inclination  on 
the  distance  of  closest  approach  is  also  evident.   This 
Figure  predicts  an  increase  in  the  distance  as  the  inclina- 
tion angle  of  the  orbit  is  increased. 

Figure  12  depicts  the  effects  of  changes  in  the  angular 
momentum,  lambda,  relative  to  the  earth  on  the  approach 
distance  to  the  moon.   As  the  angular  momentum  is  increased 
(this  corresponds  to  an  increase  in  velocity  and/or  perigee 
height)  the  distance  of  closest  approach  to  the  moon  is 
increased.   The  effect  of  the  inclination  of  the  orbit  shows 
up  in  this  Figure  as  before  to  indicate  greater  approach 
distances  for  higher  values  of  inclination. 

Figure  13  can  be  used  to  predict  two  effects  which  occur 
due  to  parameter  changes.   The  effect  of  T-,  on  the  closest 
approach  distance  is  due  to  the  variation  of  the  moon's 
position  at  launch  of  the  vehicle.   For  positive  values  of 
T,  the  moon  and  the  vehicle  interact  very  closely,  whereas 
for  negative  values  of  T1  the  moon  passes  far  from  the 
vehicle. 

The  effect  of  a  shift  on  the  distance  of  closest 
approach  may  also  be  seen  in  Figure  13.   Smaller  values  of 
the  apse  angle  allow  for  closer  interaction  with  the  moon. 
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The  value  for  planar  orbits  is  independent  of  the  apse  angle 
and  therefore  remains  the  same. 

Figures  16,  17,  and  18  show  the  effect  of  these  same 
parameters  on  the  distance  of  closest  approach  for  orbits 
of  only  slight  inclinations.   As  can  be  seen  from  these 
figures  the  effects  of  changes  in  rho,  lambda,  and  T,  are 
the  same  as  before.   As  would  be  expected  changes  in  the 
apse  angle  and  small  changes  in  the  inclination  angle  from 
the  planar  case  cause  no  significant  change  in  the  distance 
of  closest  approach  as  predicted  in  the  planar  case. 

Figures  iH   and  15  show  the  three-dimensional  geometry 
of  the  earth  orbit  and  the  moon  passage.   The  inclination 
angle,  i,  discussed  is  the  angle  between  the  earth-moon 
plane  and  the  orbit  plane.   The  apse  angle,  w,  is  the  angle 
from  the  node  line,  ft,  to  the  apse  axis.   Similar  parameters 
are  shown  for  the  moon  passage  orbit. 
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V.   CONCLUSIONS  AND  RECOMMENDATIONS 

The  method  of  uniformly  valid  asymptotic  expansions 
appears  to  have  merit  as  an  initial  analytic  tool.   The 
ability  of  this  solution  to  predict  values  of  parameters 
of  earth-moon  trajectories  with  such  a  small  expenditure  of 
computer  time  well  justifies  use  of  these  expansion  solu- 
tions.  A  quick  parametric  variational  study  could  predict 
"ball-park"  initial  conditions  which  would  greatly  reduce 
the  time  of  total  solution  of  a  trajectory  problem. 

Observing  the  results  of  the  analysis  of  this  paper  it 
appears  that  one  might  conclude  that  the  closest  approach 
possible  occurs  when  the  launch  energy  is  as  high  as  possible 
The  angular  momentum  relative  to  the  earth  must  be  low  also 
for  this  closest  approach  to  occur.   The  moon's  position 
which  is  fixed  by  T1   and  co  must  be  such  that  it  is  just 
ahead  of  the  particle  as  they  come  near  each  other. 

The  only  other  parameter  which  must  be  considered  is  the 
launch  inclination.   The  minimum  value  of  closest  approach 
is  indicated  to  occur  for  planar  trajectories. 

From  this  discussion  one  then  concludes  that  minimum 
distance  of  closest  approach  occurs  for  a  planar  launch  at 
the  highest  value  of  energy  and  lowest  value  of  angular 
momentum  possible.   These  are  attained  for  near  earth 
launch  with  high  perigee  velocities. 


M 


It  is  recommended  that  further  study  should  be  conducted 
in  the  vicinity  of  the  minimum  energy  solution  region.   This 
is  the  region  of  the  previous  mentioned  instability  and  is 
very  important  in  the  study  of  low  energy  transfers  where 
most  of  todays  travel  takes  place. 

It  would  also  be  interesting  to  release  some  of  the 
restrictions  on  the  analysis  in  order  to  apply  the  solution 
to  interplanetary  problems. 

Kevorkian  [8]  has  extended  the  analysis  to  higher  order 
terms  and  a  comparison  if  this  result  with  those  higher 
order  solutions  would  also  be  of  interest. 

Additional  work  is  necessary  before  the  present  analysis 
can  be  applied  to  a  full  spectrum  of  trajectory  problems. 
The  present  work  however  indicates  that  useful  information 
may  be  gained  from  this  form  of  solution  for  the  restricted 
three-body  problem  and  that  this  solution  has  merit  as  a 
potential  analytical  tool  for  initial  design  studies. 
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FIGURE    !.    AMES,  VAN  DYKE  BOUNDARY   LAYER  PROBLEM. 
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FIGURE  2.   FLOW  CHART   OF  MATCHING   PRINCIPLE. 


47 


Departure 
Asymptote 


Moons  Orbit 


FIGURE  3.  EXAMPLE   EARTH- MOON   HYPERBOLIC  ENCOUNTER 
TRAJECTORY. 
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FIGURE  4,     THREE -BODY  GEOMETRY 
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FIGURE  5  ,   EARTH-MOON  COORDINATE    SYSTEMS   FOR  THE 
RESTRICTED    3- BODY  PROBLEM. 
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FIGURE  6.    ANGULAR    MOMENTUM   CORRECTION. 
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FIGURE  7.  ANGULAR    MOMENTUM  RELATIVE   TO  THE  MOON. 
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FIGURE  8.     DISTANCES  OF  CLOSEST  ENCOUNTER  WITH  MOON. 
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FIGURE  9.  PERIGEE  HEIGHT   VS.  LAMBDA  FOR  VALUES 
OF  RHO. 
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FIGURE  10.   PERIGEE  VELOCITY    VS.  LAMBDA  FOR  VALUES 
OF  RHO. 
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FIGURE  II.    INCUNATION   VS.  CLOSEST   APPROACH   FOR 

VARYING  ENERGY,  ALL  OTHER    FARAMETERS 
CONSTANT 
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FIGURE  12 .  INCLINATION  VS.  CLOSEST  APPROACH  FOR 
VARYING  ANGULAR  MOMENTUM,  ALL  OTHER 
PARAMETERS    CONSTANT. 
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FIGURE   14,    EARTH   ORBIT    3-D  GEOMETRY 
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FI&RE   15.   MOON   PASSAGE    3-D  GEOMETRY 
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FIGURE  16.    APSE   ROTATION  ANGLE   VS.  CLOSEST  APPROACH 
.     FOR   3   VALUES    OF  i  WITH   p   AND  T,  CONSTANT 
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FIGURE  17.    APSE    ROTATION  ANGLE  VS.  CLOSEST 
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FIGURE  18.  APSE  ROTATION  ANGLE  VS.  CLOSEST 
APPROACH  FOR  3  VALUES  OF  i  WITH 
p   AND    X  CONSTANT. 
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LIST  OF  EQUATIONS  USED  IN  THE  IBM  SOLUTION 

The  flowchart  of  Figures  17  and  18  shows  how  the  logic 
of  the  computer  solution  is  carried  out.   The  following 
equations  are  those  used  in  the  solution  on  the  computer. 


INPUT 


j  o^b  _  Vp^   1  ''r 


A  -  -T ***' 7"- 

I  \AA.  -    _L  1 


£j  r  Apse      Lvt.A-rto+^> 

EVALUATION  OF  TERMS  FOR  X  =  1 

Cap-to  =  To   =  toot.^ 
<2   =  3   ^  Ci-hu^V^ 


en 


EVALUATION   OF   TERMS    IN    x-y    (MOON)    SPACE 

XLBA*  s    J    =-\l  +>l6l+OzSiK,H)/z  4ULV-0  -u(To-t-T^-VA.°)]AosI 
A*A£  5  6-    -      I  /zrtv 

EP5BAC-E    6    -      [/+    ^fcl^''^ 

284*    =£     -     S/AJ-'(-£  Cl-hU^/j 

EVALUATION  OF  VELOCITY  COMPONENTS 

Ul  =  .a    =  u 

VI  B     V,       ^~\ 

WU   W,     ^  o 

LUIS     0„     r^^)^^^)^-'!^^^6^'7^^1 

EVALUATION    OF   EARTH    ORBIT   PROPERTIES    AFTER   ENCOUNTER 

HEff  =  he.,  -  "fc  4  Vu  -  '/2. 
tftltr  <ge„  -   -CVn+^V„     4  I 

F&H  =    f>eit   =    dv/u    4-OU.i 
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CLOSEST  APPROACH 


COMPUTE  PARAMETER  CHECKS  OF  LAGERSTROM  AND  KEVORKIAN 

TAW   =  ucr-i)   =  UCCFH  -CC-i'i) 

?,4/2/Z  =  J  +  UT{      —  AV^ 

EQUATIONS  USED  IN  THE  SUBROUTINES  TO  OBTAIN  THE  MOTION 
INTEGRALS  .,  „  I 

Too*    ~-      l-WlWtfk-f**-)^^^*]]-*? 

TooPi   =  -took)  r  \j^l%U  -/»*•*)! 
F^       *    U)  .  C,»CWrt-T.1-*  _  G0S  &>.«  -To] 

J  ',+   ^-2/fc  Cos  OcfcO  -To]   S     *" 
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